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ABSTRACT 

Spectral indices are useful tools for quantifying the strengths of features in moderate-resolution spec- 
tra and relating them to intrinsic stellar parameters. This paper focuses on the 4300A CH G-band, 
a classic example of a feature interpreted through use of spectral indices. G-band index definitions, 
as applied to globular clusters of different metallicity, abound in the literature, and transformations 
between the various systems, or comparisons between different authors' work, are difficult and not 
always useful. We present a method for formulating an optimized G-band index, using a large grid of 
synthetic spectra. To make our new index a reliable measure of carbon abundance, we minimize its 
dependence on [N/Fe] and simultaneously maximize its sensitivity to [C/Fe]. We present a definition 
for the new index SziCH), along with estimates of the errors inherent in using it for [C/Fe] determi- 
nation, and conclude that it is valid for use with spectra of bright globular cluster red giants over a 
large range in [Fe/H], [C/Fe], and [N/Fe]. 

Subject headings: stars: abundances, stars: red giants, techniques: spectroscopic 



1. INTRODUCTION 

Spectroscopic indices are a measure of the strength of 
an absorption or emission feature, usually reported as 
the magnitude difference between the integrated flux in 
the region of the feature in question (the "science band") 
and one or two nearby continuum regions (the " compar- 
ison bands"). Indices and narrow-band photometry have 
a long history of usefulness in terms of physical infor- 
mation returned per time spent observing. They can be 
used to study a wide range of stellar properties: Ca II 
H&K surveys quickly identify extremely metal-poor stars 
(e.g., Beers et al 1985, 1999), the Mg I index is an indi- 
cator of surface gravity in cool stars (e.g., Morrison et al. 
2000), and Balmer indices are diagnostics of stellar age 
and metallicity when seen in galaxy spectra (e.g., Kauff- 
mann et al. 2003). Because they require far less observ- 
ing time than high-resolution spectroscopy, the telescope- 
limited astronomer frequently uses these low-resolution 
methods in a first pass to identify interesting objects to 
be followed up later in more detail. 

Narrow-band colors are often used to calculate phys- 
ical quantities: an empirical calibration can be estab- 
lished between the observed quantities and the intrinsic 
properties in a well-studied sample, and that calibration 
can then be applied to derive the intrinsic properties of 
a much larger set of stars, as long as the stars in the 
larger set are similar enough in temperature, gravity, and 
composition t o the well-studied calibration sample. The 
lAlonso et al.l (|1999f l photometric temperature calibration 
is a well-known example of this process. The observ- 
ables can also be used directly as proxies; for example, 



DDO and Stromgren colors, which are known to corre- 
late with temperature, metallicity, and surface gravity, 
can be used for dwarf/giant separation, or to do rough 
sorting by metallicity (e.g., Grundahl et al. 2002). The 
so-called Lick indices can be measured from integrated 
spectra of galaxies, and provide information about unre- 
solved stellar populations (e.g., Worthey 2004). 

Molecular bandstrength indices measured from low- or 
moderate-resolution spectroscopy are often the end re- 
sult themselves: comparing CN versus CH bandstrengths 
of globular cluster red giants, one finds that CH band- 
strength falls with rising luminosity (|Carbon et al.|[l982T ) 
(the same is true in the halo field: see Gratton et al. 
2000), while at a given luminosity the CH- weaker (and 
therefore lower carbon) stars tend to have the stronger 
CN bands, and therefore must be nitrogen-enhanced. 
These trends imply that red giant atmospheres must be 
experiencing ongoing CN(0)-cycle processing (e.g., Mc- 
Clure & Norris 1977 , Shetrone et al. 1999). Much re- 
search into the behavior of C and N elemental abun- 
dances in globular cluster red giants (e.g., Briley & Co- 
hen 2001, Norris et al. 1984, Norris & Zinn 1977, Zinn 
1973) has relied heavily on the use of the 4300 A CH 
bandstrength for [C/Fe] determination. This is true of 
Population II field giants as well (e.g., Rossi et al. 2005). 
However, there has been little uniformity in the CH in- 
dices used by various researchers. Some G-band indices 
that have been used to date are S(CH) (jMartell et all 
2008): 
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which was defined for use with the low-metallicity glob- 
ular cluster M53 ([Fe/H]=-1.84); 
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s C h (|Brilev fc Smithlll993h : 
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which was defined for use with the moderate-metallicity 
globul ar cluster M13 ([Fe/ H1=-1.54); 
m C H (|Trefzger et al.lll983t ): 
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which was defined for use with red giants in metal-poor 
glob ular clusters M92 and M15 ([Fe/H] ~ -2.3); 
GP (iRossi et al.ll200l : 
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defined for a study of carbon-enhanced metal-poor field 

stars; 

and Cff(4300) (lHarbeck et alJl2003h : 



CH(4300) = -2.5 log 
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which was defined for use with main-sequence stars of 
the globular cluster 47 Tucanae. In the above equations, 
I\ is the flux per A ngstrom, recorded in data such as 
iMartell et all (|2008l ) in units of ADU per pixel, and F v 
is the flux per unit f requency, calib r ated to true energy 
units in the study of iTrefzger et alj (|1983l ). 

2. DEFINING A NEW G-BAND INDEX 

We are motivated to define a G band index that is 
valid across a wide range of metallicity for use in stud- 
ies of carbon and nitrogen abundances in bright globular 
cluster red giants. Since this region of the spectrum has 
broad CN features as well as broad CH features, minimal 
dependence on nitrogen abundance is an important cri- 
terion for our index definition, in addition to predictable, 
monotonic dependence on overall stellar metallicity and 
[C/Fe]. This is doubly important in light of the work we 
intend to use this index for, converting CH and 3883A 
CN bandstrengths to carbon and nitrogen abundances 
to study the rate of ongoing CN(0)-cycle processing in 
globular cluster red giants. 

In the present analysis, we use a large grid of synthetic 
spectra generated wit h MARCS model atmospheres 
(|Gustafsson et al.lll975l ) and the SSG spectrum-synthesis 
code (Bell et al. 1994 and references therein). These 
are the same synthe tic spectra used fo r index-based 
abund ance analyses in lMartell et al.l (|2008D . ICohen et all 
(2002), and other papers by the current authors. Our 
particular grid represents bright red giants with My = 
— 1.5, twelve values of [Fe/H] ranging from —2.31 to 
-0.83, [C/Fe] varying from -1.4 to +0.4 in steps of 0.2 
dex, and [N/Fe] ranging from —0.6 to +2.0, also in steps 
of 0.2 dex. Other model parameters, identical for each 
model, are [0/Fe]= +0.2, [Ca/Fe]= +0.3, 12 C/ 13 C = 10, 
and i>turb = 2 km/s. Spectra were smoothed to a resolu- 
tion of 5.4A and a pixel spacing of 1.8 A to better match 
the data we intend to study with this index. Figure 1 




Fig. 1. — The G-band region of two synthetic spectra from our set 
with [Fe/H]= -1.41, [C/Fe]= +0.4, and [N/Fe]= -0.6 (thin line) 
and +2.0 (heavy line). The vertical dotted lines mar k the science 
band of the low-metallicity G-band index S(CH) UMartell et al.1 
120081) . The broad nitrogen-sensitive feature to the blue of the G 
band is the 4215A CN band. Spectra are normalized at 4319. 2A, 
just redward of the G band, to emphasize nitrogen-sensitive fea- 
tures. 
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Fig. 2. — Synthetic spectra demonstrating the full range of [N/Fe] 
(upper panel) and [C/Fe] (lower panel) available in our model grid. 
Spectra are normalized just redward of the G band, as in Figure 1. 
Nitrogen sensitivity is almost entirely confined to the 3883A and 
4215A CN bands, while carbon sensitivity is nearly universal. 

shows the relevant region of two of the synthetic spectra, 
with [Fe/H]= -1.41, [C/Fe]= +0.4, and [N/Fe]= -0.6 
(thin line) and +2.0 (heavy line), normalized at 4319. 2A 
to emphasize relative differences. The G b and is marked 
by v ertical dotted lines at the S{CH) (|Martell et al.l 
|2008f ) science band, and the broad nitrogen-sensitive fea- 
ture blueward of the G band is the 4215A CN band. The 
upper panel of Figure 2 shows synthetic spectra with typ- 
ical globular cluster values for [Fc/H] (—1.41 dex) and 
[C/Fe] (-0.4 dex), and the full range of [N/Fe] available 
in our model grid. The values of [N/Fe] range from —0.6 
(weakest absorption) to +2.0 (strongest absorption) in 
steps of 0.2 dex. The lower panel of Figure 2 shows an 
analogous set of spectra, with [N/Fe] fixed at +0.6 and 
[C/Fe] varying from —1.4 to +0.4 dex, the full range rep- 
resented in our models, also in steps of 0.2 dex. There is 
very little wavelength range that is unaffected by carbon 
or nitrogen in blue spectra of stars of this type. 



3 



2.5 ■ 

2.0 : 



[Fe/H] = -2.31 



• • • • • 



2.5 
2.0 



[Fe/H] = -2.31 



1.0 - 

o.5 ; 
o.o ; 
— -o.5 ; 

\ 2.5 ; 

" 2.0 ; 
1 -5 ; 
1 .0 - 
0.5 - 
0.0 - 
-0.5 ': 



• • • 

• • • 

• • • 



-1.41 
••• < 



• • • 



• • • • • • 



• 

«•••••••• • 

«•••••••• • 

I 

«•••••••• I 

*•••••••• I 

«••••••• i 



0.9 

S(CH) 



0.6 



0.7 



FlG . 3. — Low-metallicity G-band index S{CH) iJMartell et alj 
2008) versus [N/Fe] for four example metallicities spanning our 
model grid. [C/Fc] ranges from —1.4 (leftmost column in each 
panel) to +0.4 (rightmost column in each panel) in steps of 0.2 
dex. While S(CH) is fairly unaffected by [N/Fe] at low metallicity, 
there is a marked decrease in S(CH) at high values of [N/Fe] as 
[Fe/H] rises. This effect makes S(CH) an unreliable indicator of 
[C/Fe] for stars in the high abundance range. 

Figure 3 shows the CH index S{CH) (jMartell et all 
|2008i) versus nitrogen abundance for the full range of 
[C/Fe] and [N/Fe] available in our grid of synthetic spec- 
tra. In each panel, [C/Fe] increases from —1.4 (leftmost 
column) to +0.4 (rightmost column) in steps of 0.2 dex. 
The [Fe/H] metallicity varies from —2.31 (top left panel) 
to —0.83 (lower right), and as the overall metallicity in- 
creases, the resulting index-abundance grid becomes less 
rectangular. At the highest metallicity, high nitrogen 
abundance causes a large decrease in the S(CH) index 
for a fixed carbon abundance, making it difficult to deter- 
mine [C/Fe] from S(CH) without first knowing [N/Fe]. 
This effect can be attributed to the presence of CN ab- 
sorption features in the blue comparison bandpass of the 
S(CH) index. Equiv alent plots for the index CH (4300) 
(jHarbeck et al.l 120031 ). which was used for a study of 
main-sequence stars in 47 Tuc, are shown in Figure 4. 
The decrease in the measured index value seen in Figure 
3 has been overcome(by placement of the blue compari- 
son bandpass redward of the A4215 CN bandhead), and 
slightly overcompensated for, producing a more regular 
grid of points for all but the largest values of [Fe/H], 
[C/Fe] and [N/Fe]. However, the horizontal range in 
these plots is smaller than in Figure 3, meaning that 
there is a smaller change in Cff (4300) for a given change 
in [C/Fe]. Ideally, our new index will combine the good 
characteristics of these two figures, with strong depen- 
dence on [C/Fe] and independence from [N/Fe]. 

3. OPTIMIZATION 

We use the general function minimizer TNMIN 1 to 
define a new index Sz(CH) in a way that simul- 
taneously maximizes its sensitivity to carbon abun- 
dance and minimizes the influence of the nearby CN 
features. The quantity that TNMIN uses to deter- 
mine the nitrogen sensitivity of an index is ^^^/r^ > 
the mean of the absolute value of the change in 
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FlG . 4. — Main-sequence G-band index Off (4300) JHarbeck et al.1 
2003) versus [N/Fe] for the same four example metallicities and 
[C/Fe] options shown in Figure 3. The effect of [N/Fe] on measured 
bandstrcngth is smaller than in Figure 3, but the sensitivity to 
[C/Fe] is also lower. 

S2(CH) per unit dex change in [N/Fe]. This pa- 
rameter is calculated between consecutive [N/Fe] val- 
ues in the subset of model grid points that have max- 
imal [C/Fe] (e.g., 5 3(S'2(Ci/)|[c/Fc]=+0.4 : [N/Fo]=+2.0 - 

S , 2 (Cif)|[ C /Fc]=+o.4 : [N/Fc]=+i.8))- As can be seen in Fig- 
ure 3, the models with the largest [C/Fe] abundance 
(the rightmost column of points) are diverted the fur- 
thest from a rectangular grid by high [N/Fe] and [Fe/H] 
abundances. Minimizing the mean derivative between 
consecutive points in that column will force it to be as 
vertical as possible, which will mean that the definition 
for S2(CH) will be as nitrogen-independent as possible. 

A similar quantity a[c/Fc] ^ s our measure of the sen- 
sitivity of S-z(CH) to carbon abundance. It defined to 
be the difference between S^Ci?) measured at maxi- 
mal [C/Fe] and minimal [N/Fe], and S2 (CH) measured 
at minimal [C/Fe] and minimal [N/Fe], divided by the 
range of [C/Fe] available in our model grid. In Figure 

3, for example, ^c/FV|' > * s measured as the difference 
in S2{CH) between the rightmost and leftmost points of 
the lowest row, divided by 1.8 dex. A large value for 
A^cff wiU correspond to a S 2 (CH)-[N /Fe] plot cov- 
ering a large horizontal range, and will reduce the error 
inherent in using S^Cif) to determine [C/Fe]. 

We assume a priori that minimizing the effect of [N/Fe] 
on CH bandstrcngth will mean keeping the comparison 
and science bands from overlapping with the CN fea- 
tures, primarily the 4215A band, while maximizing the 
response to carbon will mean restricting the science band 
to the G band. We approached the two CH bandstrength 
criteria independently, first fixing the location of the sci- 
ence band and allowing TNMIN to find the best loca- 
tions for the comparison bands by minimizing nitrogen 
dependence, and then fixing the comparison bands and 
having TNMIN return the optimal science band loca- 
tion by maximizing carbon dependence. This process 
is then repeated, with the results of the previous opti- 
mization used as inputs for the next, until the result is 
self-consistent. The problem is solved independently for 
each of the twelve [Fe/H] values in the synthetic spec- 
trum grid, and we find that the result is not strongly 
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Fig. 5. — Nitrogen sensitivity measure ^^fj^/Kp versus central 

wavelength for a 30- A- wide blue comparison band, calculated as- 
suming a science band from 4297A to 4317A and a red comparison 
band from 4330 A to 4375 A, for the typical globular cluster metal- 
licity [Fe/H]= —1.41. Shaded regions show the two blue compar- 
ison band locations preferred by the TNMIN routine. Synthetic 
spectra with [Fe/H]= -1.41, [C/Fe]= +0.4, and [N/Fe]= -0.6 
(thin line) and +2.0 (heavy line) are also plotted (uppermost two 
curves) to show the spectral behavior inside the possible compari- 
son bands. 



dependent on overall stellar metallicity. 

There are two virtually equivalent options for the lo- 
cation of the blue comparison band returned by the TN- 
MIN program, 4040A to 4050A and 4207A to 4238A. 
Since the index we are trying to define will be used on 
data with a signal to noise ratio per pixel between 40 
and 100, we would like for the bands to be wide, so that 
pixel-to-pixel random errors will be averaged out. We 
would also like for the bands to be closely-spaced, so 
that errors in flux calibration will be minimally impor- 
tant, making the resulting bandstrength measurement 
more easily intercomparable with data taken on other 
dates, with other instruments and by other observers. 
By both of these criteria, placing the blue sideband at 
4207A to 4238A will produce a better index. To fur- 
ther explore the dependence of ^TnT^t on tri e Pl ace ~ 
ment of the blue sideband, we calculated ^^/p^ for 

all twelve [Fe/H] values available in our model grid, using 
151 definitions of S2{CH) with a red comparison band 
at 4330A to 4375 A, the S{CH) science band (4280A to 
4320A), and a 30- A- wide blue comparison band centered 
at wavelengths ranging from 3950. 5A to 4250.5A. This 
produced twelve ^^n/r^ versus ^center curves, which 
are averaged together to create the curve shown in Fig- 
ure 5, which also includes TNMIN's preferred locations 
for the blue comparison band (shaded regions) as well as 
synthetic spectra with [Fe/H]= —1.4, [C/Fe]= +0.4, and 
[N/Fe]= -0.6 (thin line) and +2.0 (heavy line). The 

two relative minima in the ^TnTr!] curve correspond 
relatively well with the TNMIN results, and with expec- 
tations: both are at the edges of the 4215A CN band. 

Figure 6 shows a subset of our models in more de- 
tail: the four lower curves are ^n/fV]^ calculated for a 

30-A-wide band centered at a wavelength varying from 
4170A to 4270A, for the same four metallicities used 
in Figure 3. The ideal location, the central wavelength 



Fig. 6. — More detailed version of Figure 5, with 



\AS 2 (CH)\ 



A[N/Fc] 

versus central wavelength for a 30-A-wide comparison band, calcu- 
lated for four sample metallicities spanning the range of our model 
grid. The location of the minimum in each curve (which produces 
the least nitrogen-sensitive G-band index) is dependent on [Fe/H]. 

where 'ArM^r^n is minimal, is dependent on [Fe/H], with 



A[N/Fc] 

higher-metallicity curves preferring 
that overlaps more with the 4215/ 



tion, the minimum value of /Ont^i 



a comparison band 
CN band. In addi- 
is also dependent 

on metallicity: the lower minima of the lower-metallicity 
curves mean decreases with decreasing metallicity, mean- 
ing that, as seen in Figures 3 and 4, 5*2 (CH) is less de- 
pendent on [N/Fe] at low metallicity. Interestingly, op- 
timal comparison band locations at all metallicities in- 
clude a small region of the 4215A CN bandhead in the 
blue comparison band, and all of the ^^N/iff curves 
level off at non-minimal values as soon as the blue com- 
parison band moves entirely out of the 4215A absorption 
band. This is unexpected given that we are explicitly 
designing S2(CH) to be independent of nitrogen abun- 
dance. The implication is that either the science band 
or the red comparison band contains a region of spec- 
trum that is sensitive to [N/Fe], and it is necessary to 
have a slightly nitrogen-sensitive blue comparison band 
to compensate. To choose a location for the blue compar- 
ison band that will work equally well for all metallicities, 
we once again average the ^^/r^ curves calculated 
for all twelve [Fe/H] values together. That mean curve 
is minimized when the blue comparison band runs from 
4211. 5A to 4241.5A. This is close to, but not exactly 
the same as, the result of our recursive TNMIN pro- 
cess. The difference comes from the use of metallicity- 
averaged ^MTpfr^ curves in this subsequent refinement 
of the TNMIN result. 

The location of the red comparison band returned by 
TNMIN is vir t ually the s ame as in th e defi nitions of 
lHarbeck et all (|2003h and iMartell et all (|2008l ): 4330A 
to 4375A. Since there is very little [C/Fe]-independent 
spectrum redward of the G band (see Figure 2), we adopt 
this result directly as the optimal location for the red 
comparison band instead of trying to improve on it. 

Fixing the comparison bands in their optimal loca- 
tions and allowing the science band location to vary, we 
find that TNMIN strongly prefers a narrow science band. 
This is reasonable; the strongest part of the CH feature 
has the strongest dependence on [C/Fe]. However, it is 
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Fig. 7. — Metallicity-averagcd a[C/Fc ] (middle panel) and 
^^N/Fcj^ (lower panel) curves versus wavelength Xbi ue of the blue 
end of the science bandpass. The top panel plots synthetic spectra 
for [Fc/H]= -1.41, [N/Fe]= +2.0, and [C/Fc]= -1.4 (thin line) 
and +0.4 (heavy line) overplottcd to show the location of the G 
band. Each individual curve is calculated by fixing the red edge 
of the S2(CH) science band at a value between 4312A and 4321A, 
then varying the blue edge from 4230A to 4310A. Comparison 
bands are held fixed at their optimal locations (4212A to 4242 A 
and 4330A to 4375A). As described in the text, the red edge of the 
science band is set based on the ^ ^a [ n/FcP ^ curves, and the blue 



edge location is based on the 
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impractical for the science band to be too narrow; as with 
the comparison bands, we want to have a broad enough 
wavelength coverage that noise in the spectrum will av- 
erage out to zero. To better understand the optimal 
placement of the S 2 (CH) science band, we calculated 

AS (CH) 

the carbon- and nitrogen-sensitivity measures A [c/ Fc ] 

and ^j^/F^ for a wide variety of possible S2(CH) 
index definitions. In each case, the comparison bands 
were fixed at their optimal locations, as discussed above. 
For ten values of X re d (the position of the red edge of 
the science band) between 4312A and 432lA, we calcu- 
lated ^jffiffi and '^ffiffi 1 for 161 possible values of 
Xbiue (the position of the left edge of the science band) 
and all twelve [Fe/H] values available in our model grid. 
We then collapsed the twelve independent [Fe/H] curves 
corresponding to each value of X re d into ten metallicity- 
averaged curves. We base our determination of the opti- 
mal values for X re d and \u U e on these mean curves. 

The ten metallicity-averaged sensitivity curves are 
shown in the middle and lower panels of Figure 7 as a 
function of Xbiue- In the upper panel synthetic spec- 
tra are shown as a function of A for [Fc/H]= —1.41, 
[N/Fe]= +2.0, and [C/Fe]= -1.4 (thin line) and +0.4 
(heavy line). Generally, lower values of X re d result in 

larger values for both ^[c^Fc] anc ^ ^n/Fc^ ' though 
the ^n^T^rn curves with larger X re d values show an 



A[N/Fo] 

upturn for values of Xu% 



above 4300A. The rise in 



^Alc^Fc] as tne k lue e< ig e °f the science band approaches 
the strongest part of the G band is a clear visual rep- 
resentation of TNMIN's preference for a narrow science 
band for S 2 (CH). 



Fig. 8. — In analogy to Figures 3 and 4, this plot shows S2(C'H) 
versus [N/Fe] for four sample metallicities and ten [C/Fc] values 
spanning our model grid. The decrements in index value seen at 
high [Fc/H] and [N/Fe] in Figure 3 have been reduced at moder- 
ate [Fc/H] and replaced by a slight increase at high [Fe/H], and 

the horizontal range ^gyp^" 1 is larger at all metallicities than in 

Figures 3 and 4. 

Our criteria for the science band, in order of impor- 
tance, are a low value for ^^jp^-, a high value for 

^c/Fc^ ' anc ^ an overan bandwidth larger than 15 A. We 
select X re d by requiring that ^j^/p^ be lower than 

0.01 in the range 4290A< Xuue < 4300A, then choos- 
ing the smallest value of X re d that qualifies to allow for 
the largest possible ^c/Fc^ ■ This results in a value 
for X re d of 4318A. To choose Xbiue, we simply find the 
highest point on the ^gyp^ curve corresponding to 

Xred = 4317A, which gives a value for Xbi ue of 4297A. 
This result satisfies our bandstrength criterion. We then 
performed a check for self-consistency, using this final 
science band location as the input for placement of the 
blue comparison band. We found that the resulting com- 
parison band was nearly identical to our original result 
(4212 A to 4242 A). Repeating the process of fixing the 
comparison bands to choose a science band, we find that 
it also shifts slightly. We adopt the range 4297A to 4317A 
as the optimal location for the science band of S2(CH), 
producing a final bandstrength definition of 
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Figure 8 is analogous to Figures 3 and 4, and shows 
S 2 {CH) versus [N/Fe] for the full [C/Fe] and [N/Fe] 
range available in our grid of synthetic spectra, at four 
example metallicities. The points in this plot, in par- 
ticular those at high [C/Fe] and high [N/Fe], are ar- 
ranged much more uniformly and perpendicularly than 
those shown in Figures 3 and 4. The horizontal range 
covered is also larger than in Figures 3 and 4, meaning 
that S 2 (CH) is more sensitive to [C/Fe] than S(CH) or 
CH (4300). As a result, the transformation from mea- 
sured S 2 (CH) to [C/Fe] is more straightforward and less 
error-prone. 

The lower right panel of Figure 8, showing S 2 (CH) 
versus [N/Fe] for synthetic spectra with [Fe/H]= —0.83, 
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Fig. 9. — Method-based measurement error A[C/Fe] ver- 
sus [Fe/H], [C/Fe], and [N/Fc] for 10 4 randomly-chosen ([Fe/H], 
[C/Fe], [N/Fe]) abundance points within the range of our model 
grid. The error A [C/Fe] is a complex function of abundance, 
but is never prohibitively large: 99.7% of points with typical 
globular cluster red giant abundances (—2.31 < [Fe/H] < —0.83, 
-1.0 < [C/Fe] < 0.0, 0.0 < [N/Fe] < +1.0) fall in the range 
-0.03 < A[C/Fe] < +0.03. 

does show deviation from rectangularity at large [C/Fe] 
and [N/Fe] values. Fortunately, in the stars this in- 
dex is designed for, this abundance combination is 
rarely, if ever, observed to exist. In bright red gi- 
ants in globular clusters, large nitrogen abundances are 
strongly correlated with low carbon abundances (in gen- 
eral, [C+N+O/Fe] is roughly constant in these stars, e.g., 
Denissenkov et al. 1998 and Sneden et al. 1997), though 
the mechanism for producing these particular abundance 
ratios is a subject of ongoing study. As a result, this un- 
avoidable imperfection in our index definition should not 
present a large problem for studies of bright red giants 
in globular clusters. 

4. APPLICATIONS AND LIMITATIONS OF THIS INDEX 
DEFINITION 

In order to characterize the error introduced by us- 
ing S^iCH) with an assumed nitrogen abundance to de- 
termine [C/Fe], we performed a series of Monte Carlo 
tests, in each trial selecting 10 4 random ([Fe/H], [N/Fe], 
[C/Fe]) abundance sets within the parameter space oc- 
cupied by our synthetic spectra and then interpolating 
within the model grid to find the corresponding "ob- 
served" values of ^ (Cff). We used a p rocess similar to 
the one described in lMartell et al.l ()2008l ) to convert each 
randomly-chosen value of ^(Cff) to [C/Fe], first inter- 
polating within the model grid to match the randomly- 
chosen [Fe/H], then assuming a [N/Fe] value of +0.6 
and interpolating linearly within that subgrid to find the 
[C/Fe] value needed to produce the "observed" S^Ci?). 
The difference between the calculated and original input 
[C/Fe] values is the error A [C/Fe]. We also performed 
this test with the extra step of adding Poisson-distributed 
1%, 2%, 5% and 10% noise to the synthetic spectrum cor- 
responding to the randomly-chosen abundance set and 
redetermining the value of ^(Cff), to understand how 
much additional error in the calculated [C/Fe] might be 
introduced by noise in the data. 

Figure 9 shows A [C/Fe] versus [Fe/H], [N/Fe], and 
[C/Fe] for the no- noise test, and there are two things 
immediately obvious from these figures. First, that the 




-1.5 -1.0 



[Fe/H] 



2.0 -1.5 



Fig. 10.— Each panel shows A [C/Fe] versus [Fe/H] for 10 4 
randomly-chosen ([Fe/H], [C/Fe], [N/Fc]) abundance points within 
our model grid. We added Poisson noise to the spectra correspond- 
ing to those abundances, then calculated [C/Fe] and A [C/Fe] based 
on the value of S2 (CH) measured from the noisy spectra. The noise 
ranges from 1% (upper left) to 10% (lower right). While the shape 
of the A [C/Fe] distribution remains the same in all panels, the 
magnitude is strongly dependent on the amount of noise added to 
the spectrum. The standard deviation 0A[C/Fe] roughly doubles in 
each successive panel. 

behavior of A [C/Fe] is complex and not simply related 
to any one parameter. Second, S2{C'H) is fairly reliably 
converted to [C/Fe]: the absolute value of A [C/Fe] is 
never larger than 0.11 dex, and the vast majority of er- 
rors are quite small. The mean A [C/Fe] value is —0.0049, 
and 80% of all trials were in the range -0.03 < A [C/Fe] 

< 0.03. 

Those numbers change to —0.024 and 58% in trials 
with [N/Fe]> 1.0, to -0.0031 and 72% in trials with 
[Fe/H]> -1.4, and to 0.0046 and 78% in trials with 
[C/Fe] > 0.0. In the region of abundance space typically 
occupied by globular cluster red giants (—2.31 < [Fe/H] 

< -0.83, -1.0 < [C/Fe] < 0.0, 0.0 < [N/Fe] < 1.0 ), the 
mean A [C/Fe] value is 0.0020, and 99.7% of trials have 
-0.03 < A [C/Fe] < 0.03. 

Figure 10 shows A [C/Fe] as a function of [Fe/H] for 
one particular Monte Carlo realization of each Poisson 
noise-added test. In tests with 1% and 2% noise added 
to the spectra, the S2(CH)-to-[G/Fe] conversion pro- 
cess is still quite reliable. Mean values for A [C/Fe] are 
-0.0089 and 0.0098, respectively, and 66% and 49% of 
points, respectively, fall in the range —0.03 < A [C/Fe] 

< 0.03. Considering only the region of abundance space 
likely to be occupied by globular cluster red giants, those 
numbers change to —0.0008 and 84% (for 1% error) and 
-0.0004 and 60% (for 2% error). The tests with 5% 
and 10% error are less encouraging: while the shape of 
the A [C/Fe] distribution remains the same in all pan- 
els of Figure 10, the standard deviation <7a[c/Fc] roughly 
doubles with each increase in noise. Again restricting 
the abundance range to the region where globular clus- 
ter red giants are mostly likely to be found, the mean 
value for A [C/Fe] and percentage of stars falling in the 
range -0.03 < A [C/Fe] < 0.03 are -0.0067 and 27%, re- 
spectively, for the 5%-error case and —0.0067 and 14%, 
for the 10%-error case. This implies that [C/Fe] abun- 
dances derived from St(CH) (or any other C-band in- 
dex) should only be trusted when the original spectra 
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Fig. 11. — A comparison of [C/Fel values d e termi ned using 
S 2 (CH), and [C/Fe] values from ISmitn~B rilcv (2006), for four 
giants in M13 observed by the authors as part of an ongoing study. 
While there is a fairly predictable relationship between the two sets 
of measurements, there is clearly an offset between the zeropoints 
of the two carbon abundance scales. 

have decently high signal to noise ratios: A [C/Fe] is un- 
acceptably large in the 10% error case, and is marginal 
in stars with typical globular cluster abundances in the 
5% error case. 

It is important to test our new index on data, and 
not just noise-injected synthetic spectra, to demonstrate 
that it does return reliable [C/Fe] abundance values. 
Figure 11 shows a comparison of [C/Fe] values derived 
from 52 (CH) with [C/Fe] values taken f rom the litera- 
ture compilation of lSmith fc Brilevl ((2006) (originally re- 
ported in Suntzeff 1981), for four M13 giants observed by 
the present authors as part of an ongoing study of carbon 
depletion in globular cluster red giants. The spectra were 
taken with the Kast double spectrograph on the 3-meter 
Shane telescope at Lick Observatory, and have moder- 
ate resolution (R~ 1000). The synthetic spectra were 
smoothed to match the data resolution and pixel spac- 
ing. Carbon abundances were derived from measured 
S2{CH) bandstrengths through a n index-matching pro - 
cess similar to the one described in lMartell et al.l (2008). 

Figure 11 sho ws a clear linear relationship between the 
ISuntzeffl (|1981h [C/Fe] abundances and the [C/Fe] val- 
ues based on S2(CH) for the four stars in common be- 
tween the samples, though the abundance scales defined 
by the two different abundance-determination methods 
have different zeropoints. The magnitude of the offset be- 
tween the two abundance scales i s quite similar to that re- 
ported inlSmith & Brilevl (|2006D . roughly —0.35 dex. In 
ISmith fc Brilevl (|2006l) the au t hors f ind that abundance- 
scale offset between ISuntzefj (|1981h [C/Fe] abundances 
and those r eported in studies usin g MARCS model at- 
mospheres (jGustafsson et al.l I1975T) . Our carbon abun- 
dances based on 5*2 (CH) can therefore be readily com- 
pared with other MARCS-based abundance determina- 
tio ns, and r e liably transformed onto the abundance scale 
of ISuntzefil (|1981h and compared with other measure- 
ments made on that scale. 

Figure 12 shows method-based errors for carbon abun- 
dances as a function of [Fe/H], analogous to the leftmost 
panel of Figure 9, but expanded to show the same re- 
sult calculated using the indices S(CH) and sch- Both 
of these indices were defined for studies of one particu- 
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Fig. 12.— Each panel shows A[C/Fe] versus [Fe/H] for 10 4 
randomly-chosen ([Fe/H], [C/Fe], [N/Fc]) abundance points within 
our model grid. In the leftmost panel, A [C/Fe] was calculated from 
measurements of S2(CH); in the center, from S(CH), and on the 
right, sch- The vertical scale in all three panels is identical. 

lar globular cluster, S(CH) for the low-metallicity clus- 
ter M53 and sch for the moderate-metallicity cluster 
M13, and therefore return small errors at and below those 
clusters' respective metallicities. The leftmost panel in 
Figure 12 is identical to the leftmost panel in Figure 
9, and the vertical range is the same in all three pan- 
els. The mean value in the central panel of Figure 12, 
the difference between input and measured carbon abun- 
dances calculated using S(CH), is -0.035, and 69% of 
trials fall in the range -0.03 < A[C/Fe]i < 0.03. In 
the rightmost panel, which shows the difference between 
input and measured carbon abundances calculated from 
sch, the mean value is 0.007, and 86% of trials have 
-0.03 < A[C/Fe] 2 < 0.03. The advantage of the S 2 (CH) 
index over these other two is that it has been designed 
to be fairly insensitive to [Fe/H] and [N/Fe], resulting in 
a smaller overall range of A[C/Fe] errors across a wide 
range of [Fe/H], [N/Fe] and [C/Fe]. 

5. CONCLUSIONS 

It is our hope that this new G-band index can provide 
a standard tool for low-resolution spectroscopic studies 
of bright red giants in globular clusters. It is simultane- 
ously more sensitive to [C/Fe] and less sensitive to [N/Fe] 
than similar indices used in the literature to date, and 
would allow for useful comparison and combination of in- 
dependent researchers' data in multiple globular clusters. 
The optimization method used here could be generally 
applicable to any spectral region where there are a few 
main parameters controlling the strength of an absorp- 
tion or emission feature. Our characterizations of the 
dependence of S 2 (CH) on [Fe/H], [C/Fe], and [N/Fe], 
as well as the error associated with converting 5*2 (CH) 
to [C/Fe], apply over a wide abundance range, but they 
are only strictly true for stars with My = —1.5, since 
the model atmospheres and synthetic spectra were only 
calculated for that particular value. While our current 
model calibration of 5*2 (CH) should lead to a reliable 
measure of [C/Fe] for a moderate range in My, perhaps 
— 1.0 > My > —2.0, it should not be applied to stars out- 
side this magnitude range without recalibration. Since 
the index was designed specifically for bright red giants 
of Population II metallicities, we do not recommend us- 
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ing it to study [C/Fe] in stars belonging to significantly 
different populations. 
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